Lie- Semigroup Structures for Reachability and Control of Open Quantum Systems: 
Kossakowski-Lindblad Generators Form Lie Wedge to Markovian Channels 



G. Dirr/'B U. Helmke,^ I. Kurniawan/ and T. Schulte-Herbruggen^'0 

^Institute of Mathematics, Wurzburg University, Am Hubland, D-97074 Wiirzhurg, Germany 
^Department of Chemistry, Technical University Munich, 
Lichtenbergstrasse 4, D-85747 Garching, Germany 
(Dated: 22nd February 2009) 

In view of controlling finite dimensional open quantum systems, we provide a unified Lie-semigroup 
framework describing the structure of completely positive trace-preserving maps. It allows (i) to 
identify the Kossakowski-Lindblad generators as the Lie wedge of a subsemigroup, (ii) to link 
properties of Lie semigroups such as divisibility with Markov properties of quantum channels, 
and (iii) to characterise reachable sets and controllability in open systems. We elucidate when 
time-optimal controls derived for the analogous closed system already give good fidelities in open 
systems and when a more detailed knowledge of the open system (e.g., in terms of the parameters of 
its Kossakowski-Lindblad master equation) is actually required for state-of-the-art optimal-control 
algorithms. As an outlook, we sketch the structure of a new, potentially more efficient numerical 
approach explicitly making use of the corresponding Lie wedge. 
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I. INTRODUCTION 

Understanding and manipulating open quantum sys- 
tems and quantum channels is an important challenge 
for exploiting quantum effects in future technology [l|]. 

Protecting quantum systems against relaxation is 
therefore tantamount to using coherent superpositions 
as a resource. To this end, decoherence-free subspaces 
have been applied bang-bang controls have been 
used for decoupling the system from dissipative interac- 
tion with the environment, while a quantum Zeno ap- 
proach Q may be taken to projectively keep the system 
within the desired subspace 0. Very recently, the op- 
posite approach has been taken by solely expoiting re- 
laxative processes for state preparation [a, l7||. It is an 
extreme case of engineering quantum dynamics in open 
systems [l|jwhere targeting fix points has lately become 
of interest [3| . 

In either case, for exploiting the power of system and 
control theory, first the quantum systems has to be char- 
acterised, e.g., by input-output relations in the sense of 
quantum process tomography. Deciding whether the dy- 
namics of the quantum system thus specified allows for a 
Markovian description to good approximation (maybe up 
to a certain level of noise) has recently been addressed 
[lol . [ill . This is of crucial interest, since a Marko- 
vian equation of motion paves the way to applying the 
power Lie-theoretic methods [l3 . [l^ from geometric and 
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bilinear control theory. Moreover, it comes with the well- 
established frameworks of c omp letely positive semigroups 
and Kraus representations (TEElia! 18, 19, 20, 21, 22]. 

On the other hand, the specific Lie-semigroup aspects 
of open quantum systems clearly have not been elabo- 
rated on in the pioneering period 1971-76 of completely 
positive semigroups [l3, E3, [H, [13, [IS] , mainly since ma- 
jor progress in the understanding of Lie semigroups was 
made in the decade 1989-99 [H, [IB, [1^, [13, [2i,l2§| • While 
relations of Lie semigroups and classical control theory 
were soon established, e.g., in [H, [13, [sil, [13, [13, [H, [sl, 
[3^ . only recently the use of Lie-semigroup terms in the 
control of open quantum systems was initiated [stI . [131, 
where in [37| the elaborations were confined to single two- 
level systems. However, we see a great potential in ex- 
ploiting the algebraic structure of Lie-semigroup theory 
for practical problems of reachability and control of open 
quantum systems. 

Its importance becomes evident, because among the 
generic tools needed for the current advances in quan- 
tum technology (for a survey see, e.g., [l|]), quantum 
control plays a major role. From formal description of 
quantum optimal control [13| the theoretical aspects of 
existence of optima soon matured into numerical algo- 
rithms solvin g pr actical problems of steering quantum 
dynamics (iol. l4l[. [H, [i3|- Their key concern is to find 
optima of some quality function like the quantum gate fi- 
delity under realistic conditions and, moreover, construc- 
tive ways of achieving those optima given the constraints 
of an accessible experimental setting. For a recent in- 
troduction, see [i^l. However, realistic implementations 
in open quantum systems are mostly beyond analytical 
tractability. Hence numerical methods are often indis- 
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pensible, where gradient-like algorithms are the most ba- 
sic, but robust tools. Thus they proved applicable to 
a broad array of problems including optimal control of 
closed quantum systems lisl. [45| and computing entan- 
glement measures [i^ . [ZtI. l48l| . For mathematical details 
on gradient systems as numerical tools for constrained 
optimisation, we refer to (iol. Isol. IsH . 

Generalising these well-established gradient tech- 
niques, in our previous work (i^l, we have exploited the 
geometry of Riemannian manifolds related to Lie groups, 
their subgroups, and homogeneous spaces in a common 
framework for setting up gradient flows in closed quan- 
tum systems. There we addressed (a) abstract optimisa- 
tion tasks on smooth state-space manifolds and (b) dy- 
namic optimal control tasks in the specific time scales of 
an experimental setting. Here, we will see that the corre- 
sponding abstract optimisation tasks for open quantum 
systems are much more involved, while the dynamic op- 
timal control tasks remain in principle the same. From 
a mathematical point of view, this difficulty results from 
the fact that the evolution of a controlled open quantum 
system is no longer described by a semigroup of unitary 
propagators, i.e. by a semigroup contained in a compact 
Lie group. 

Thus, we extend the Lie-theoretic approach in [Hi to fi- 
nite dimensional open quantum systems and discuss their 
dynamics in terms of Lie semigroups. In particular, we 
characterise the Lie properties (the Lie wedge) of Marko- 
vian quantum channels from the viewpoint of divisibil- 
ity and local divisibility in semigroups. — On a general 
scale and with regard to practical applications of quan- 
tum control, knowing about the Lie-semigroup structure 
of the dynamic system is shown to be highly advanta- 
geous: analysing its tangent cones (Lie wedges) allows 
for addressing problems of reachability, accessibility, con- 
trollability and actual control in a unified frame providing 
powerful Lie algebraic terms. 



Starting Point 

To begin with, we briefly indicate how the theory eluci- 
dated in previous work [40| can be extended to reachable 
sets of non necessarily controllable systems. In partic- 
ular, we concentrate on the structure of reachable sets 
and obstacles arising from it. Moreover, pertinent appli- 
cations to open relaxative quantum dynamical systems 
are elaborated — proving the relevance of the semigroup 
setting in physics. 

The starting point in ji^ was a smooth state-space 
manifold M or a controllable dynamical system on AI, 
i.e. a control system whose reachable sets Reach(Xo) sat- 
isfy Reach(Xo) = M for all Xq £ M. For a right invari- 
ant system (0]) the state space of which is given by a 
connected Lie group G, controllability is equivalent to 
the fact that the entire group G can be reached from the 



unity 1, i.e. 

G = Reach(l) (J Reach(l, T), (1) 

T>0 

where Reach(l, T) denotes the reachability set in time 
T > 0, i.e. the set of all states to where the systems can 
be steered from 1 S G in time T, cf. Eqn.®. In general, 
however, we cannot expect Eqn.|(l|) to hold. Neverthe- 
less, the reachability sets Reach(l, Ti) and Reach(l, T2) 
of right invariant systems obey the following multiplica- 
tive structure 

Rcach(l, Ti) • Reach(l, T2) = Rcach(l, Ti + T2). 

Thus Rcach(l) is a subsemigroup of G, see Sec lII Dl — 
Now, we will give a basic survey on subsemigroups and 
some of their applications in quantum control. 

II. FUNDAMENTALS OF LIE 
SUBSEMIGROUPS AND REACHABLE SETS 

A. Lie Subsemigroups 

For the following basic deflnitions and results on Lie 
subsemigroups we refer to [^i . [13, [H, [H, \E^. However, 
the reader should be aware of the fact that the termi- 
nology in this area is sometimes inconsistent. Here, we 
primarily adopt the notions used in [23|. For further 
reading we also recommend (36| . 

A subsemigroup of a (matrix) Lie group G with Lie 
algebra g is a subset S C G which contains the unity 1 
and is closed under multipHcation, i.e. S • S C S. The 
largest subgroup contained in S is denoted by E(S) := 
S n S~^. The tangent cone of S is deflned as 

L(S) := {7(0) I 7(0) = 1, 7(t) e S, i > 0} C 0, 

where 7 : [0, 00) G denotes any smooth curve con- 
tained in S. In order to relate subsemigroups to their 
tangent cones, we need some further terminology from 
convex analysis. A closed convex cone tn of a flnite di- 
mensional real vector space is called a wedge. 

Moreover, a wedge ro in a Lie algebra g is termed a Lie 
semialgebra if the wedge ir is locally compatible with the 
Baker-Campbell-Hausdorff (bch) multiplication X*Y := 
X + Y + ^[X,Y] + . . . , deflned via the BCH series. More 
precisely, there has to be an open bch neighbourhood 
S C g of such that ro is locally invariant under *, i.e. 

(ro ns) * (tr n S) C m. (2) 

For a thorough treatment of the BCH multiplication and 
Lie semialgebras see [13, H^l ■ 

The edge of tt) denoted by E(rD) is the largest subspace 
contained in tn, i.e. one has E(tTi) := inn (—in). Finally, a 
wedge tn of a finite dimensional real (matrix) Lie algebra 
2 is called a Lie wedge if it is invariant under the group 
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of inner automorphisms Inn(tTi) := (cxp(adE(rD))). More 
precisely, 

for all g £ E(tt3). Here and in the sequel, we denote by 
(M) and {M)s the group and, respectively, semigroup 
generated by the subset Af C G. 

Remark II. 1. While every Lie semialgebra is also a Lie 
wedge, the converse does in general not hold, as will be of 
importance in the paragraph on divisibility in Sec. Ill CI 

Now, the fundamental properties of the tangent cone 
L(S) can be summarised as follows. 

Lemma II. 1. Let S be a closed subsemigroup of a Lie 
group G with Lie algebra Q and let Vo <Z Q be any Lie 
wedge. Then the following statements are satisfied. 

(a) The edge of Vo, E(tTi), carries the structure of a Lie 
subalgebra of q. 

(b) The tangent cone L(S) coincides with 

L(S) = {.g € I cxp(tg) G S for all t > 0}. (3) 

Ln particular, L(S) is a Lie wedge of g which is 
^d-E(s)-invariant, i.e. GxoG^^ = vo for all G e 
E(S). 

(c) The edge of L(S) fulfills the equality E(L(S)) = 
L(E(S)). 

Proof. 

(a) Note that e**^^" (/i) g E(tr) for alH G K and g, h G 
E(rD). Hence 

^e*^'^«W|,^„ = ad3/iGE(ir) 

for all g,h <E E(rD), thus E(rD) is a Lie subalgebra. 

(b) The proof of Eqn. JSl) is rather technical and there- 
fore we refer to [2J|, Proposition IV. 1.21. Once 
Eqn. ^ is established, one has 

L(S) = p|t-iexp-i(S) 

t>o 

and thus the continuity of the exponential map im- 
plies that L(S) is closed. To see that L(S) is a 
wedge we have to show: (i) AfL(S) = L(S) for all 
M G M+ and (ii) L(S) + L(S) C L(S). Property 
(i) is obvious; property (ii) follows by the Trotter 
product formula 

e*(»+'') = lim (e*^/"e*'*/")". 

n—^oo 

Finally, let g G E(L(S)) and h G L(S), then 

for all t > 0. Thus eSheS = e^^9(/i) G L(S). The 
same argument applies to G G E(S). 



(c) Let g G E(L(S)). Then e*^ G S for all t G R. Thus 
6*9 G E(S) and hence g G L(E(S)). Therefore, 
we have shown E(L(S)) C L(E(S)). The converse, 
L(E(S)) C E(L(S)), holds by definition. 

For more details, see Proposition 1.14 in [IS]. ■ 

For closed subsemigroups. Lemma IH.ll provides the 
justification to call the tangent cone L(S) Lie- or Lie- 
Loewner wedge of S. 

Unfortunately, the 'local-global-correspondence' be- 
tween Lie wedges and (closed) connected subsemigroups 
is not as simple as the correspondence between Lie sub- 
algebras and Lie subgroups. On the one hand, there are 
Lie wedges Vo such that 'the' corresponding subsemigroup 
S is not unique, i.e. the equality tn = L(S) holds for more 
than one subsemigroup S. On the other hand, there are 
Lie wedges ir which do not act as Lie wedge of any sub- 
semigroup, i.e. Vo = L(S) fails for each subsemigroup S, 
cf. 

Another subtlety in the theory of semigroups arises 
from the fact that there may exist elements in S that are 
arbitrarily close to the unity but do not belong to any 
one-parameter semigroup completely contained in S (a 
standard example being a certain subsemigroup of the 
Heisenberg group [H, Ea])- This somewhat striking fea- 
ture arises whenever the bch multiplication leads outside 
the Lie wedge L(S). It does not occur as soon as L(S) 
also carries the structure of a Lie semialgebra, cf. The- 
orem III. 21 below. In general, however, the exponential 
map of a zero-neighbourhood in L(S) need not give a 
1-neighbourhood in the semigroup. 

Meanwhile, the following terminology is well- 
established (29l . [55| : a set i? G G is called exponential 
if to each element T ^ E there exists a Lie algebra 
element g & Q such that cxp{g) = T and cxp{tg) G E 
for all t G [0,1]. Now, let S be a closed subsemi- 
group of a Lie group G with Lie wedge L(S) and let 
(expL(S))5 := {e^i •••e^" | gt G L(S), n G N} be the 
subsemigroup generated by expL(S) C G. Then 

(i) S is called Lie subsemigroup if it is characterised 
by the equality S = (cxpL(S))g,; 

(ii) S is called weakly exponential if exp L(S) is dense 
in S, i.e., if S = cxpL(S); 

(iii) S is called exponential if the set S is exponential in 
the above sense, i.e., if S = expL(S); 

(iv) S is called locally exponential if there exists a 
1-neighbourhood basis with respect to S consisting 
of exponential subsets. 

The inclusions expL(S) C expL(S) C (expL(S))g are 
obvious. A Lie wedge vo is said to be global in G if there 
exists a Lie subsemigroup S C G so that L(S) = tti, i.e. 
S = (exp(tr))s. 
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Remark II. 2. For the sake of completeness note that 
the term Lie subsemigroup is closely related (with subtle 
distinctions) to the notions of (completely or strictly) in- 
finitesimally generated subsernigroups, which will not be 
pursued here any further, cf. |24| - 

B. The Reductive and the Compact Case 

Based on the classical Cartan decomposition of reduc- 
tive Lie groups [HBl, we reformulate a known result on 
the existence of global Lie wedges — a setting which does 
arise in open quantum systems, cf. Theorem IIILSl and 
Corollary IIILGl below. We do so by stating a convenient 
version of a more gene ral result, cf. Theorem V.4.57 and 
Remark V.4.60 in [2j|, streamlined here in view of prac- 
tical application. 

Theorem II. 1. Let G be a closed connected (matrix) 
Lie group which is stable under the conjugate transpose 
inverse, i.e. which is invariant under the involution Q : 
X (X^^Y . Let g = t p be the decomposition of its 
Lie algebra into +1 and —1 eigenspaces of the involution 
De(l) ^-.e-.X^ -Xt. Then 

(a) the map p x K ^ G, (p, K) ^ exp{p)K with K := 
(expt) is a diffeomorphism onto G; 

(b) the set S := exp(c) ■ K. is a Lie subsemigroup with 
L(S) — c (B i, provided c C p is a closed pointed 
cone, i.e. E(c) = {0}. 

Proof. Combining Proposition 7.14 in [HHl with the 
proof of Theorem V.4.57 in [131, the result follows. ■ 

Fortunately, the somewhat intricate general scenario 
just outlined simplifies dramatically when considering 
compact Lie subsemigroups. 

Proposition II. 1. (24l . l36l ]. Let S be a compact .sub- 
semigroup of a Lie group G. Then S itself is a compact 
Lie subgroup ofG. 

C. Divisibility and Local Divisibility in Semigroups 

Here, we briefly summarise some results on divisibility 
in semigroups that will be useful in Section IIIICI when 
relating them to recent findings by Wolf et al. on the 
divisibility of quantum channels. 

For semigroups, there is the following well-established 
notion of divisibility [13, US]: a subset of Z) C G is termed 
divisible, if each element T € D has roots of any order 
in D, i.e. to any r S N there is an element S <E D with 
= T . Similarly, a semigroup S is called locally divisi- 
ble, if there is a 1-neighbourhood basis in S consisting of 
divisible subsets. 

For linking global and local notions of divisibility with 
exponential semigroups. Lie semialgebras play a crucial 
role. Here we start with some basic results before sketch- 
ing what became known as 'i/ie divisibility problem\ For 



details see the literature given in Further Notes and Ref- 
erences below. 

Proposition II. 2. |24!| A closed subsemigroup S of a 
connected Lie group G is divisible if and only if it is 
exponential, i.e. cxpL(S) = S. 

Proof. If S = cxpL(S), then S is trivially divisible. 
The converse is already more technical to show and we 
refer to Theorem V.6.5 in [131 . B 

Theorem II. 2. For a closed semigroup S the following 
assertions are equivalent: 

(a) the Lie wedge of S is a Lie semialgebra; 

(b) S is locally exponential; 

(c) S is locally divisible. 

Proof. For the equivalence (a) 4=> (c) see Corol- 
lary 3.18 as well as [2J| Propositions IV. 1.31-32 and Re- 
mark IV. 1.14. While the implication (6) =^ (c) is trivial, 

(a) =^ (b) follows by [121 Proposition 3.17(a). For a sim- 
ilar result on Lie semigroups see also [1^ Theorem III. 9 
and III.21. ■ 

The difficulty to go beyond the straightforward results 
just mentioned made the following closely related ques- 
tions notorious as Hhe divisibility problem^ [il. [29I. \El\ : 

(i) Is the Lie wedge ir = L(S) of a closed divisible 
i.e. exponential semigroup also a Lie semialgebra? 

(ii) When does (global) divisibility imply local divisi- 
bility? 

These problems were open for several years until settled 
in the sterling monography by Hofmann and Ruppert in 
1997 [2§|, where all Lie groups and subsemigroups with 
surjective exponential map are classified. — For studying 
local divisibility in the connected component of the unity 
in more detail (and in view of follow-up work), some of 
its main results can be summerised as follows. 

Theorem II. 3. |29l | Let G be a connected Lie group 
containing a weakly exponential subsemigroup S with Lie 
wedge Yo = L(S). 7/ S is closed and has non-empty inte- 
rior in G and its only normal subgroup is 1 S G, then 

(a) S is divisible (exponential), i.e., expL(S) = S; 

(b) its Lie wedge in = L(S) is a Lie semialgebra; thus 

(c) S is also locally divisible (locally exponential). 

Proof. For (a) see Theorem 7.3.1 and Scholium 7.3.2 
in [2§| (p 132) lifting Eggert's work [2^ on Lie semialge- 
bras to reduced weakly exponential subsemigroups thus 
leading to Theorem 8.2.14 in [2§| (p 152); assertion (b) is 
Theorem 8.2.1(v) in ^ (p 145); finally (c) follows from 

(b) by virtue of Theorem III. 21 above. ■ 
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Further Notes and References. — A (somewhat jerry- 
built) primer on divisible semigroups including an ac- 
count of earlier results and problems can be found in 
[stI . while the current status is documented in [2^. A 
broad overview on historical aspects of a Lie theory of 
semigroups is given in j5l,[5§|. Ultimately, readers inter- 
ested in links to Hilbert's Fifth Problem and topological 
semigroups are referred to (60| . 

D. Reachable Sets 

Let (S) be a right invariant control system 

X = AuX, A„ e 0, ueU C R" (4) 

on a connected Lie group G with Lie algebra g and let 
5 C 2 denote its system Lie algebra, i.e. s := {A^ \ u S 
h{)uc is by definiton the Lie subalgebra generated by Au, 
u £ U. The reachable set Reach(Aro) of (S) is defined 
as the set of all AT S G that can be reached from ATo 
by an admissible control function u{t). More precisely, 
let Xu[t) denote the unique solution of Eqn. (0]) which 
corresponds to the control u{t). Then 

Reach(A:o) := |J RGach(Xo,T) 

T>0 

with 

Reach(Xo, T) {Xu{T) G G | T > 0, u{t) e U]. (5) 

Moreover, (E) is called accessible, if Reach(Xo) has non- 
empty interior in G for all Xq S G, and controllable, if 
Reach(Aro) = G for all A'o € G. For more details on 
the control theoretic terminology and setting we refer to, 
e.g., [l3,[3ll,[6l|. Now, in the following series of results the 
relation between reachable sets of right invariant control 
systems and subsemigroups will be clarified. 

Theorem II. 4. (l4l . l36| |. Let (E) be a right invariant 
control system on G given by Eqn. Then the follow- 
ing statements are equivalent: 

(a) The system (E) is accessible. 

(b) The reachable sei Reach(l) is a subsemigroup of G 
with non-empty interior. 

(c) The entire Lie algebra q of G is generated by 
Au, u eU, i.e. s = g. 

Theorem II. 5. (36l |. Lei (E) be a right invariant control 
system on a connected Lie group G given by Eqn. ^ 
and assume that (E) is accessible, i.e. s = q. Then the 
following statements are satisfied: 

(a) The closure of the reachable set Rcach(l) is a Lie 
subsemigroup ofG, i.e. 

S= (expL(S))s 



where S := Rcach(l). Moreover, 
int S = int (Reach(l)), 

and 



S = Reacho(l), (6) 

where Reacho(l) denotes the reachable set of the so- 
called extended system, i.e. the system where A^ is 
allowed to range over the entire Lie wedge L(S). 

(b) The sei L(S) is the largest subset of g satisfying ^ 
and, moreover, it is the smallest Lie wedge which 
is global in G and contains A^, u eU. 

In control theory, due to the characterisation given in 
part (b) of Theorem IILS^ the Lie wedge L(S) is usually 
known as the Lie saturate of Au, u €U, see, e.g., 
Ig^ . Conversely, one has the following result. 

Theorem II. 6. |36l |. Let G be a connected Lie group 
and let S be a Lie subsemigroup ofG. Then, there exists 
a right-invariant control system (E) on G with control 
set {Au I u e Z^} eg such that 



S := Rcach(l). 
In particular, one may choose {Au \ u G U} = L(S). 

Finally, we summarise some well-known necessary and 
sufficient controllability conditions for right invariant 
control systems. While the first criterion is rather diffi- 
cult to check, as the computation of the global Lie wedge 
corresponding to a given control set Au is in general an 
unsolved problem, the second one provides a simple al- 
gebraic test for compact Lie groups, cf. Proposition IIL II 

Corollary II. 1. Let (E) be an accessible right invariant 
control system on a connected Lie group G, i.e. s — g. 
Then the following statements are equivalent: 

(a) The system (E) is controllable. 

(b) The Lie wedge o/Reach(l) is all of q. 

Proof. The implication (a) =^ (b) is trivial; the con- 
verse (b) (a) follows from Theorem in.4r b) and The- 
orem Hit a), cf. (lii. ■ 

Corollary II.2. [13, Let (E) be a right invari- 

ant control system on a connected compact Lie group G. 
Then controllability of (E) is equivalent to accessibility, 
i.e. to 5 ~ Q. 

Remark II. 3. If the assumption s = g in Theorem III. 51 
and Corollarv lII.il is not fulfilled, the above results, how- 
ever, still remain valid when restricting to the unique Lie 
group Go := (exps). 
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III. DEVELOPMENTS IN VIEW OF 
APPLICATIONS TO QUANTUM CONTROL 

A. Reachable Sets of Closed Quantum Systems 

An application of Corollary III.2I to closed finite- 
dimensional quantum systems, e.g., n spin-i qubit sys- 
tems with possibly non-connected spin-spin interaction 
graph yields an explicit characterisation of their reach- 
able sets. The same result based on a sketchy controlla- 
bility argument can be found in j63|. 

Theorem III.l. Assume that the spin-spin interaction 
graph, which corresponds to the controlled n spin-\ sys- 
tem 

n 

U=-i(Hd+ ^kHk^o)u (7) 

k=\ 

with Hd := J2k<i •^kicrk,zO-Lz and Hk^a ■= CFk,a, a G 
{x,y}, decomposes into r connected components with nj 
vertices in the j-th component. Then, the reachable set 
Reach(l2n) of Eqn. ^ is given (up to renumbering) by 
the Kronecker product SU{2"-^ ) S'[/(2"-). 

Proof. Suppose that the spin-^ particles of the system 
are numbered such that the first component of the graph 
contains the vertices 1, . . . , ni, the second one the vertices 
ni -\- \, . . . ,ni -\- n2 and so on. Thus n = -\- ■ ■ ■ -\- n^. 
Then, it is straightforward to show that the system Lie 
algebra is equal to the Lie algebra of Go := 5f7(2"i) ® 
• • •(»5'C/(2"-) cf. [sil. Therefore, we can consider Eqn. ^ 
as a control system on Gq. Since Go is a closed subgroup 
of SUiT'^), it is compact and thus Corollary [IL2] applied 
to Gq yields the desired result. ■ 

Henceforth read iV 2" for n spin-i qubits. — Note 
that the same line of argument as above applies to the 
modified control term discussed in (63| . 

B. Open Quantum Systems and Completely 
Positive Semigroups 

In open relaxative quantum systems jH, [H, [sl, [IB, [fi^l 
however, the situation is different because relaxation 
translates into 'contraction'. Thus the dynamics on den- 
sity operators is no longer described by the action of a 
compact unitary Lie group as before. 

Moreover, we use the following short-hand for the total 
Hamiltonian 

Hu:= Hd + ^UjHj, (8) 

j 

where uj and Hj denote possibly time dependent control 
amplitudes and time- independent control Hamiltonians, 
respectively. Now, we consider a finite dimensional con- 
trolled Master equation of motion 

p=-i adjf„ (p) - r(p) ^-Cu{p), ueUc K" (9) 



on the set of density operators 

po5i(7V) := {p G 0l(iV,C) I p = p^ p > 0, trp = 1} 

modelling a finite dimensional relaxative quantum sys- 
tem. Here, ad^^ denotes the adjoint operator, 
i.e. adij^(p) := [Hu,p], and —T represents the infinitesi- 
mal generator of a semigroup {exp(— tF) \ t > 0} oi linear 
trace- and positivity-preserving (super-)operators [lOd^ . 
Clearly, £„ and thus Eqn. |[9]) extend to the vector space 
of all Hermitian matrices 

l)cx{N) := {H e Ql{N,C) \ H = i^t}. 

Now it makes sense to ask for the self-adjointness of 
r with respect to the Hilbert-Schmidt inner product 
{Hi,H2) := iv[HiH2) on [)er(iV). Unfortunately, T need 
not be self-adjoint, yet it is self-adjoint, e.g., if it can be 
written in double-commutator form, cf. Eqn. (|22|) . 

Moreover, since the fiow of Eqn. ([9]) is trace preserving, 
the image of F is contained in the space of all traceless 
Hermitian matrices 

f)ero(iV) := {H e fll(A^,C) | H = H\tv H = 0}. 

Therefore, the restriction of rlfij^^j-^-) yields an opera- 
tor from \)zXq{N) to itself and thus Eqn. ^ can also 
be regarded as an equation on ^zx.q{N). To distinguish 
these two interpretations of E qn. ^ , we call the latter 
homogeneous Master equation |lOlj . Note that the ho- 
mogeneous Master equation completely characterises the 
dynamics of the open system, once an equilibrium state 
p* of Eqn. Q is known. More precisely, if £„(p,) = 
for all u G (e.g., choose p* = -^Iw for unital equa- 
tions) the dynamics of po := p — p* is described by the 
homogeneous Master equation. Finally, we associate to 
Eqn. ((9]) a lifted Master equation 

X = -CuoX, X(0) = id (10) 

on GL{{)zx{N)) and GL{\)zXq{N)), respectively. Equation 
(fTO|) will play a key role in the subsequent subsemigroup 
approach. 

For a constant control u{t) = u, the formal solution of 
the lifted Master equation Eqn. ifTOj) is given hy Tu{t) := 
exp(— tjC„). Thus 

{Tu{t)\t>Q} (11) 

yields a one-parameter semigroup of linear operators act- 
ing on ^zx{N). Actually, the operators T„(i) form a con- 
traction semigroup of positive and trace preserving linear 
operators on \)(.x{N) in the sense that 

||r„(t)(A)||i < \\A\\, 

for all A G f)er(iV), cf. [H, Recall that the trace 
norm ||A||i of ^ G f)er(iV) is given by 

N N 
i i 
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where at and Xi denote the singular values and eigen- 
values of A, respectively. The semigroup ifTTj) is said 
to be purity- decreasing if moreover all T^it) constitute 
a contraction with respect to the norm induced by the 
Hilbert-Schmidt inner product, i.e. if 

(T„(t)(p),T„(f)(p)) < {p,p) 

holds for all p £ p0Si{N) and all t > 0. In general, T„(t) 
is not purity-decreasing. However, if F is in Kossakowski- 
Lindblad form, cf. Eqn. ifTSj) . a necessary and sufficient 
condition for being purity-decreasing is unitality of F^, 
i.e. Tl{1n) = 0, cf. [63|- Thus for a unital Kossakowski- 
Lindblad term Tl, the subsemigroup 

Pi::^{Tu{t)\t>0,ueU)s (12) 

generated by the one-parameter semigroups ifTTj) is con- 
tained in a linear contraction semigroup of a Hilbert 
space. 

Remark III.l. Let Hhe a complex Hilbert space with 
scalar product (•,•). Then the linear contraction semi- 
group of Ti. is defined by 

C{n) :={T e GL{H) \ {Tv, Tv) < {v, v) for all v G H}. 

Note that Co(7i) — the connected component of the unity 
in C(?i) — is in fact a Lie subsemigroup. This is evident 
from the polar decomposition T = PU, because PU S 
C{H) with U unitary and P = P^ positive definite holds, 
if and only if the eigenvalues of P are at most equal to 
1. Thus 

Coin) = cxp(-c) ■ u{n) , 

where c denotes the cone of all positive semidefinite el- 
ements in 0[(7Y) and C/(7Y) the corresponding unitary 
group. Similarly, one can define contraction semigroups 
for real vector spaces, cf. [27| . 

Next, we briefiy fix the fundamental notion of complete 
positivity for open quantum systems. Recall that a linear 
map Tu{t) is completely positive, if r„(t) and all its ex- 
tensions of the form Tu{t) ® Im are positivity-preserving, 
i.e. 

(r„(i)®l„)(poSi(7V-TO)) Cp05i(7V-m) 

for all m S N. Complete positivity of the Markovian 
semigroup {Tu{t) | i > 0} is required to guarantee that 
{Tu{t) \ t > 0} can be associated with a Hamiltonian 
evolution on a larger Hilbert space, cf. (23I. IgqI. [7o| . 

According to the celebrated work by Kossakowski [l^ 
and Lindblad [2o|, Eqn. ^ generates a one-parameter 
semigroup {Tu{t) | i > 0} of Hnear trace-preserving and 
completely positive operators, if and only if F^ can be 
written as 

i E ^kVkP + pVlVk - 2Vt.pVl Tl{p) (13) 

k 



with arbitrary complex matrices S qI{N, C). Thus the 
Master equation ^ then speciaHses to the Kossakowski- 
Lindhlad form 

Cu{p) := i ad^/„ ip) + ^J2 ^kVkP + pV^Vk - 2VkpVl 

k 

(14) 

Suppose we consider the complexification of \)zx{N), 
i.e. the complex vector space 

\)zx{N f - Q{{N, C) = C^^^ ^ C^' . 

By extending the linear operators adu^j^L £ Q^^^'^i^)) 
to Hu,Tl ■ C^' C^" one arrives at the superoperator 
representations 

Hy^ := 1n®Hu- HJi^In and (15) 

tl kY.^N®vx + ykV:^iN-2v:^v, , (le) 
fc=i 

where Hu,Tl G 0l(Af^,C) are iV^ x N"^ complex matri- 
ces. In particular, if F^ is self-adjoint, the correspond- 
ing matrix representation F^, G g[(A^^,C) is Hermitian. 
Moreover, note that the matrix representation F^ con- 
tains some redundancies on gl{N'^,C) since the original 
Tl operates on the real vector space [)er(-/V) which has 
obviously smaller (real) dimension than C^" . Viewed in 
this way, note that F^ is not the same as the matrix 
representation of Tl in the coherence-vector formalism. 
See [g^I for an introduction on coherence vectors in open 
systems and [7l| for a recent characterisation of positive 
semidefiniteness in terms of Casimir invariants. More ge- 
ometric features can be found in [t^I- 

Now, the previous semigroup theory allows to interpret 
the Kossakowski-Lindblad master equation in terms of a 
Lie wedge condition. We define P to be the semigroup of 
all positive, trace preserving invertible linear operators 
on [)er(iV), i.e. 

P {T G GL{\)ev{N)) \ T ■ p05^{N) C poSi(iV)}. 

and P'^P to be the closed subsemigroup of all completely 
positive ones, i.e. 

pep := £ p I r completely positive} C P. 

Then, Pq and Pg"^ denote the corresponding connected 
components of the unity. Moreover, an arbitrary linear 
trace preserving completely positive, not necessarily in- 
vertible operator on \)^x{N) is usually called a quantum 
channel. Thus in terms of quantum channels, P'^p is the 
set of all invertible quantum channels. Now, a key-result 
by Kossakowski and Lindblad can be formulated as fol- 
lows. 

Theorem III.2. (Kossakowski, Lindblad [H, [2ol |) 

The Lie wedge L(PqP) is given by the set of all linear 
operators —C of the form C :— iad/^+F^, where F^ is 
defined by Eqn. Hfl^ . 
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While the finite-dimensional version of Theorem IIII.2I 
stated above was originally proven by Gorini, Kos- 
sakowski and Sudarshan at the same time Lind- 

blad [iOl handled the explicitly infinite-dimesional case 
of a norm (uniform) continuous semigroup of completely 
positive operators acting on a VF*-algebra. (Note that 
Kossakowski-Lindblad-type equations with time depen- 
dent coefficients were analysed, e.g., by [t^ or [73|.) 

For proving Theorem IIII.2[ a former, actually infinite- 
dimensional result by Kossakowski on one-parameter 
semigroups of positive (not necessarily completely pos- 
itive) operators on trace-class operators Bi{Ti.) and 
their infinitesimal generators was recast into a finite- 
dimensional setting in [l^. Although Kossakowski and 
Lindblad exploited different methods from functional 
analysis, a crucial point in both papers and is 
the theory of dissipative semigroups on Banach spaces, 
cf. Lumer and Phillips [75| . 

Yet in the context of finite-dimensional Lie semigroups, 
the same results now show up as a consequence of a more 
general invariance theorem for convex cones: roughly 
spoken the infinitesimal generator of a one-parameter 
semigroup leaving a fixed convex cone invariant is char- 
acterised via its values at the extreme points of the cone, 
cf. Theorem 1.5.27 in jl^l. In particular, Kossakowski's 
work on one-parameter semigroups of positive oper- 
ators then turns out to be a special application of the 
afore-mentioned invariance theorem to the convex cone 
of all positive semidefinite N x TV-matrices 

pos(iV) {H e l)ex{N) \ H > 0} . 



Likewise, Theorem lin.2l can be obtained by the invari- 
ance theorem applied to the cone pos(iV^), once the 
equivalence of complete positivity of exp(— 1£) and posi- 
tivity of ex p(— t C ® Im) is established, cf. [l^. For more 
details see 17611. 



C. Lie Properties of Semigroups versus 
Markov Properties of Quantum Channels 

Recall the notation P'^p for the closed semigroup of 
all completely positive invertible maps, whose connected 
component of the unity is termed Pg"^. Having derived 
the Lie wedge of Pg^, the issue of its glohality naturally 
emerges. Since Pg^ is closed in GL(\\zx{N)), an affirma- 
tive answer to this problem is obtained by Proposition 
V.1.14 in m. 



Theorem III. 3. The semigroup 



T:=(cxp(L(PS'^))),CPSP 



(17) 



generated by L(PqP) is a Lie subsemigroup with the Lie 
wedge L(T) = L(PgP). In particular, L(Pg'^) is a global 
Lie wedge. 

Ultimately, the question arises whether Pg'' is itself a 
Lie subsemigroup in the sense of Section[lll However, the 



identity T = Pg^ one might surmise is disproven by the 
fact that there are indeed invertible quantum channels T 
with detr > that do not belong to the subgroup T, 
cf. 

For relating these references to our context, we have 
to establish some of the terminology of Holevo [t^I and 
Wolf et al. [lH: Similar to our definition in Sec- 
a quantum channel T is called (inifinitely) 



tion in CI 

divisible if for all r € N there exists a channel S such 
that T = S"^. [NB: In stochastics and quantum physics 
[13, [m, [73, uE, [li it is long established to use the 
term 'infinitely divisible', whereas in mathematical semi- 
group theory it is equally long established to simply say 
'divisible' instead (see also Section [lICp . This is why 
here we use the brackets.] In contrast, a channel is said 
to be infinitesimal divisible if for all £ > there is a se- 
quence of channels Si, S2, ■ ■ ■ , Sr such that ||S'j — id|| < e 



and n, 



=1 



T. Moreoever, a quantum channel is 



termed time (in) dependent Markovian if it is the solu- 
tion of a Master equation X = —Co X , with initial con- 
dition X{Q) = id and time (in) dependent Liouvillian ~C 
of Kossakowski-Lindblad form. Now, for our purpose the 
results in (lol. can be resumed as follows. 

Proposition III.l. (Tol. [tqI] 

(a) The set of all time independent Markovian channels 
coincides with the set of all (infinitely) divisible and 
invertible channels. 

(b) The closure of the set of all time dependent Marko- 
vian channels coincides with the closure of the set 
of all infinitesimal divisible channels. 



The proof of Proposition IIII.ll (a) is given in [73|, part 
(b) is precisely Theorem 16 of 10||. Thus in relation to 
the work of Wolf et al. Theorem IIII. 31 reads: 

Corollary III.l. The closure of the set of all time de- 
pendent Markovian channels forms the Lie subsemigroup 
T defined in lfT7|) . Its tangent space at the unity is given 
by the Lie wedge L(Pg'^) of all Kossakowski-Lindblad gen- 
erators. 

However, one also arrives at the no-go result: 

Theorem III. 4. |lOl | The semigroup Pq^ is neither (in- 
finitely) divisible nor infinitesimal divisible. In particu- 
lar, there are invertible quantum channels which are not 
infinitesimal divisible. 

For N = 2, the above assertion is rigorously proven by 
Theorem 24 in [l3|. For N > 2, the statement currently 
presupposes one may extrapolate from the numerical re- 
sults (also for N = 2) i n [Tr |. 

Now, from Theorem IIII. 41 we conclude: 

Corollary III. 2. Pg^ itself is not a Lie subsemigroup. 
Yet in particular the semigroup P'^p of all invertible quan- 
tum channels is made of three subsets, all of which also 
occur in the connected component Pg^ ; 
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(a) the set of time independent Markovian channels 
which is given by definition as the union of all one- 
parameter Lie semigroups {exp(— £t) \ t > 0} with 
—C in Kossakowski-Lindblad form; 

(b ) the closure of the set of time dependent Markovian 
channels which coincides with the Lie semigroup T 
defined by ((T7|) ; 

(c) besides, there is a set of non-Markovian channels 
(i.e. neither time independent nor time dependent 
Markovian) whose intersection with Pq^ has non- 
empty interior. 

Clearly, Markovian channels of type (a) are a special 
case of type (b) and (a) i s eve n a proper subset of (b), 
since T is not exponenti al [ l02l| . There are also quantum 
channels with dot T < jlO| , but they can only occur out- 
side the connected component Pq^ , and thus they are ob- 
viously non-Markovian. The geometry of non-Markovian 
channels seems to be well-understood in the single-qubit 
case {N = 2), yet remains to be analysed in full detail 
for larger N. 

Corollary III. 3. (a) The semigroup Pq^ is neither lo- 
cally divisible nor locally exponential. 

(b) The Lie wedge ^(Pq'') of all Kossakowski-Lindblad 
generators does not form a Lie semialgebra. 

Proof. Again, for N = 2, part (a) follows from The- 
orem 24 in [10||. For N > 2, the assertion extrapolates 
from the numerical results in [ll|. Part (b) is an imme- 
diate consequence of part (a) and Theorem III. 21 ■ 

Now, the distinction between Lie wedge and Lie- 
semialgebra structure can be exploited to separate be- 
tween time dependent Markovian quantum channels and 
time independent ones. In general, this separation is 
rather dehcate. Clearly, as soon as a time dependent 
channel T has a representation of the form T = 0^=1 
such that the Si, S2, ■ . ■ , Sr generate an exponential Lie 
semigroup, then T is actually time independent. Though 
almost a tautology, this statement is quite difficult to 
check and therefore an (infinitesimal) condition that is 
easier to verify is most desirable. The following corol- 
lary is meant as a first result in this direction — with the 
shortcoming that it applies to channels close to unity. 

Corollary III. 4. Let T be a time dependent Markovian 
channel that allows for a representation T = Ilj^^i 
with Si ~ e^^^,S2 = e^^'^,...,Sr = e^'~'' and where 
ror denotes the smallest global Lie wedge generated by 
Li, C2, ■ . ■ , Cr- Then 

(a) T boils down to a time independent Markovian 
channel, if it is suffiently close to the unity and if 
there is a representation so that the associated Lie 
wedge XOr also carries Lie- semialgebra structure; 

(b) conversely, if T is a time independent Markovian 
channel, a representation with tr,. being a Lie semi- 
algebra trivially exists. 



Proof. The result follows by the same line of argu- 
ments as Corollary nil. 51 below. ■ 

Thus in summary three elucidating results have 
emerged: (i) the set of all time dependent Markovian 
quantum channels forms a Lie subsemigroup T and (ii) its 
Lie wedge coincides with with the set of all Kossakowski- 
Lindblad operators: it is the Lie wedge to the subsemi- 
group Pg"^ of all invertible quantum maps. Moreover, 
(iii) the border from time dependent to time indepen- 
dent Markovian quantum channels is characterised by the 
existence of an associated Lie wedge that specialises to 
Lie-semialgebra structure. 

D. Effective Liouvillians 

In physical applications a frequent task amounts to 
describing the evolution of a controlled Master equation 

X = oX, ueUdMT, (18) 

cf. Eqns. (jH [9]) with in Kossakowski-Lindblad from 
Eqn. lfT4|). by an appropriate one-parameter semigroup. 
More precisely, given an admissible time dependent con- 
trol u{t) € U and a final time toff > one is inter- 
ested in an effective time-independent Liouvillian CcS 
such that the two time evolutions coincide at teff > 0, 
i.e. T„(t)(tcff) = e"*""^"". This is a natural extension 
from average Hamiltonian theory of closed systems to 
average Liouvillians of open ones 

[ia,[8il,[8i,[8i. 

Now, Lie-semigroup theory provides a useful frame- 
work to settle the question under which conditions not 
only the final point e"*""'''"''' , but also the entire trajec- 
tory {e~*'~°"\Q < t < tcs} up to the final point complies 
with the Master equation (fTSl defining the physics of the 
system. 

Corollary III. 5. Given a Master equation (fTSl) and the 

smallest global Lie wedge tn generated by the set of con- 
trols {-Cu\u eU <Z K"}, cf Theorem^L^ Then the 
following assertions are equivalent: 

(a) The Lie wedge Vo also is a Lie semialgebra. 

(b) Any solution of ifTSj) coincides at least locally, 
i.e. for sufficiently small t > with some one- 
parameter semigroup generated by an effective Li- 
ouvillian CcS G tr. 

Proo f. Follows from the fact that the Lie semigroup 
(exp tTi)g is locally exponential if and only if its Lie wedge 
is a Lie semialgebra, cf. Theorem III. 21 ■ 

Only if the effective Liouvillian is guaranteed to remain 
within the Lie wedge W associated to the controlled Mas- 
ter equation ifTSj) then it generates a one-parameter semi- 
group 

{e-t£ctf \ t >0} that can be considered 'physical' 
at all times t > 0. Otherwise, the physical validity of the 
time evolution described by the semigroup {e"*^"''' |t > 0} 
is in general Hmited to a set of discrete times (including 
t = and t = toff)- 
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E. Controllability Aspects of Open Quantum 

Systems 

Structural Preliminaries 

Studying reachable sets of open quantum systems sub- 
ject to a controlled Hamiltionian, cf. Eqn. (flO)) below, 
is intricate, as will be evident already in the following 
simple scenario: consider a Master equation in the su- 
peroperator form 

vec p = - (i ^ Hj + Tl) vcc p , 

3 

where the iHj are skew-Hermitian, while shall be Her- 
mitian. Thus they respect the standard Cartan decompo- 
sition of 0[(A^^, C) := J © p into skew-Hermitian matrices 
(6) and Hermitian matrices (p). Then the usual commu- 
tator relations C 6, [p,p] C t, [6, p] C p suggest that 
double commutators of the form 

[[i?„fi],[i/fc,fi]] 

generate new {-directions in the system Lie algebra as 
will be described below in more detail. 

For the moment note on a general scale that such con- 
trolled open systems thus fail to comply with the stan- 
dard notions of controllability: not only does this hold 
for operator controllability of the lifted system but also 
for usual controllability on the set of all density opera- 
tors, cf. (stI . [38| . Hence it is natural to ask for weaker 
controllability concepts in open systems. 

For simplicity, we confine the subsequent considera- 
tions to unital systems of Kossakowski-Lindblad form, 
i.e. rL(ljv) = 0, as their dynamics is completely de- 
scribed by the homogeneous Master equation 

p = -iadH„(p)-ri(p) = -/:„(p) (19) 

on [)erg(A^) and its lift 

X = ~Cu o X (20) 

to GL{i]zx-Q{N)). Here the controlled Hamiltionian takes 
the form of Eqn. ([8]) with Hd and Hj in su(A^) and no 
hounds on the controls Uj G K. Thus the semigroup Pe 
given by Eqn. l(20| will be regarded as a subsemigroup 
o{ GL{\)zx.q{N)) in the sequel. Alternatively, by the pre- 
viously introduced superoperator representation, we can 
think of Pe as embedded in GL{N^, C). 

If, in the absence of relaxation, the Hamiltonian system 
is fully controllable, we have 

{iHd, iHj I J = 1, ... , m)ue = 5u{N) , (21) 

or, equivalently, 

{iHd, iHj I J = 1, ... , m)ue = psu(iV) C suiN^) , 



where we envisage psu(iV) to be represented as Lie 
subalgebra of su(A^^) given by all malices of the form 
H - H^ ® 1) with iH e 5u{N). Master equations 
which satisfy Eqn. l(2T|) are expected to be generically 
accessible, i.e. their system Lie algebras generically meet 
the condition 

{i&dHa +rL,iad//^ | j = 1,2, . . . ,m}ue = 0l(f)ero(iV)) , 

cf. [3^. [tgI. [s^I . Here, the system Lie algebra of the control 
system (cf. Section IH Pp is not to be misunderstood as 
its Lie wedge, which in general is but a proper subset of 
the system Lie algebra. 

The group generated by Eqn. (|20| therefore gener- 
ically coincides with GL{1)cVq{N)). Thus already the 
coherent part of the open system's dynamics, i.e. the 
'orthogonal part' of the polar decomposition of elements 
in Ps, has to be embedded into a larger orthogonal 
(unitary) group than of the same system being closed, 
i.e. when = 0. This can easily be seen if the Master 
equation (flO)) speciaHses so that the respective matrix 
representations iHj for i adj^^. are skew-Hermitian, while 

Fl is Hermitian. For instance, this is the case in the 
simple double-commutator form 

p= -(iadH„+i^ad^J(p) . (22) 

fe 

It exemplifies the details why iterated commutators 
like [[i/j, Fi], [i/fe, Fi]] typically generate new skew- 
Hermitian directions in the system Lie algebra of 
Eqn. I|20p . This holds a forteriori if — as henceforth — 
we allow for general Kossakowski-Lindblad generators no 
longer confined to be in double-commutator form (|22|) . 
We can therefore summarise the above considerations as 
follows. 

Resume. In open quantum systems that are fully con- 
trollable for Tl = 0, one finds: 

1. Only if T L\i)cxo[N)) acts as scalar 7I and thus 
[iHj^Ti,] = for all j, the open dynamics is con- 
fined to the contraction semigroup (0, 1]-Ads(7(jv) of 
the unitary adjoint group Adsu{N)- Moreover, the 
contractive relaxative part and the coherent Hamil- 
tonian part are independent in the sense that their 
interference does not generate new directions in the 
Lie algebra. 

2. Yet in the generic case, the open systems' dynam- 
ics explore a semigroup larger than the contraction 
semigroup of the unitary part KAsu{ n) of the closed 
analogue. 

Thus for an explorative overview, the task is three-fold: 
(i) find the system Lie algebra 

Sopen := (iadff^ -t-Fi, i adi/^. )ue ; (23) 
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(ii) if Sopcn = 9Kh^^o{N)) already (as will turn out to 
be the case in most of the physical applications with 
generic relaxative parts r^), then the dynamics of 
the entire open system takes the form of a con- 
traction semigroup contained in GL(l)erQ(-/V)); the 
relaxative part interferes with the coherent Hamil- 
tonian part generating new directions in the Lie 
algebra, where the geometry of the interplay deter- 
mines the set of explored states; 

(iii) in the (physically rare) event of Sopon ^ Ql{i)erQ{N)) 
the system dynamics takes the form of a contrac- 
tion semigroup contained in a proper subgroup of 
GL(f)ero(7V)). 



Weak Hamiltionian Controllability 

As mentioned before, controllability notions for open 
systems weaker than the standard one are desirable, since 
Eqn. lfT9|) is in general non-controllable in the usual sense. 
Here, we define a unital open quantum system to be 
Hamiltonian controllable (n-controllable) if the subgroup 
{Ad(7 I U € SU{N)} is contained in the closure of the 
subsemigroup Ps, i.e. 

{Adu \ U £ SU{N)} (zPj:. 

In constrast, we will call a system to be weakly Hamil- 
tonian controllable (WH-controllable) if the subgroup 
{Adu I U E SU{N)} is contained in the closure of the 
subsemigroup E+ • Ps C Gi(f)ero(-/V)), i.e. 

{Adu \ UeSU{N)}c [l,oo)-Ps. 

So far, WH-controllability has not been studied in the lit- 
erature, although it provides a partial answer to the prob- 
lem of finding the best approximation to a target density 
operator pf by elements of the reachable set Reach(p), 
where pp itself is contained in the unitary orbit 0{p). 
For establishing a first basic result on WH-controllable 
systems, the subalgebras generated by the controls terms 

ic ■= {iHi, . . . ,iHm)uc 

and by the Hamiltionian drift plus controls terms 

td ■= (i-ffd, i-ffi, • • • , ^H,n)uc 

will play an essential role. 

Proposition III. 2. A unital open quantum system ifTOj) 
with the Hamiltonian given by Eqn. ^ is 

(a) u- controllable, if tc = su{N) and no bounds on the 
control amplitutes uj , j ~ 1, . . . ,m are imposed; 

(b) WM- controllable, if Id = su(A^) andTL\^^^^ |.^j=7l 
with 7 > 0. 



Moreover, for U G SU{N), the smallest A G R"*" such that 
Adu S APs is given 

hye<T'{u)^ jj^hereT*{U) denotes the 
optimal time to steer the lifted system given by Eqn. l(20|) 
without relaxation, i.e. for = 0, from the identity 1 to 
Adu ■ In particular, for tc = 5u{N) one has A = 1 for all 
U e SU{N). 

Proof, (a) First, suppose tc ~ su{N). Then, for 
Fi = the fact that we do not assume any bounds on 
the controls Uj € M implies that one can steer from the 
identity 1 to any Adu arbitrarily fast. Thus for F^ ^ 
a standard continuity argument from the theory of ordi- 
nary differential equations shows that one can approxi- 
mate Ad;7 up to any accuracy by elements of Ps- Thus 
H-controllability holds. 

(b) Suppose td = su{N) and Fl||^^^ = 7I. By Corol- 

larvlll.2l we obtain controllability of {Adu \ U & SU{N)} 
for Fx, — 0. Therefore, we can choose a control u{t) which 
steers the identity 1 to Adu in optimal time T*(U). Ap- 
plying the same control to the system under relaxation 
yields a trajectory which finally arrives at e"''''^ Adu- 
Thus WH-controllability holds for A = e''^*^^^ Moreover, 
by the time optimality of T*{U) it is guaranteed that 
A = e'^^'(^) is the smallest A G R+ such that Ad^ e APs 
holds. ■ 

In general, an open quantum system that is fully con- 
trollable in the absence of relaxation will not be neces- 
sarily WH-controllable when including relaxation, even 
though it may be accessible. A counterexample showing 
this fact for the simplest two-level system and simula- 
tions will be provided in [tHI- Establishing necessary and 
sufficient conditions for WH-controllability of open quan- 
tum systems is therefore an open research problem. For 
unital systems which are controllable in the absence of 
relaxation, we do expect that the 'ratio' of the Hamilto- 
nian and the relaxative drift term completely determines 
WH-controllability. — Finally we will see that additional 
assumptions ensuring the preconditions of Theorem III. II 
allow for inclusion of the global Lie wedge of Eqn. I|19p . 

Theorem III. 5. Assume that the unital Master equation 
(fT9l) with the Hamiltonian given by Eqn. ^ fulfills the 
following condition: there exists a pointed cone c in the 
set of all positive semidefinite linear operators on \)tx^{N) 
such that 

2. [c, c] C a.d^u{N) and [c, a.dsu{N)] C c - c ; 

3. Adc/ c Adcz-i c c for all U e SU(7V) . 

Then, the Lie subsemigroup Ps of Eqn. (fT9|) is contained 
in the Lie subsemigroup 

exp(-c) ■ AdsuiN) 

with Lie wedge (— c) ® adjufjv). 




Figure 1: (Colour online) Quantum state-space manifolds for open relaxative systems shown as subsets of t)cxQ(N) with scales 
corresponding to the metric induced by the Hilbert-Schmidt scalar product. The centre of the high-dimensional sphere is the 
zero-matrix, and the geometry refers to larger systems, e.g., multi-qubit systems with A'' > 4. If in the absence of relaxation, 
the system is fully controllable, the reachable set for a fixed initial state represented as density operator po takes the form of 
the entire unitary orbit Adsu(N)ipo)- It serves as a reference and is shown as closed curve in red. In the text we focus on two 
different scenarios of open systems: (a) Dynamics of weakly Hamiltonian controllable systems with the Kossakowski-Lindblad 
term acting approximately as scalar Vl — 7I are confined to the subset (marked in blue) of states evolving from po under the 
action of the contraction semigroup (0, 1] ■ KAsu{n)- The latter is depicted as grey surface of a 'funnel' intersecting the surface 
of the high-dimensional sphere in the unitary orbit. Towards the origin, i.e., at long times, the reachable set of WH-controUable 
systems typically wraps the entire surface (dark blue portion), (b) In the generic case when [Fl, H^] / («/ = d; 1, 2, . . . , m), 
the dynamics with initial state po evolves within the volume shown in blue. New directions due to the interplay of coherent 
Hamiltonian evolution and relaxation make the dynamics explore a much larger state space than resulting from the simple 
contraction semigroup (0, 1] ■ Adgjjjjv), i.e. the surface in part (a) or even the volume contained in its interior. The intersection 
(green portion) of the volume Reach(po) with the surface of the sphere consists of the set of all states reachable from po in zero 
time or without relaxative loss. This may often collapse to the single point po or its local unitary orbit [ssl. Is^. 



Proof. By Theorem lILST b). it is sufficient to verify 
that exp(— c) • Ad^j/^^) is a Lie subsemigroup with Lie 
wedge (— c) ©adsu(7v)- This will be achieved by applying 
Theorem IILII To this end, we define g := 6 © p with 
t := &d^xx{N) and p := (c — c) + (c — c)^. Note that 
the set c — c consisting of all differences within the cone 
c coincides with the vector space spanned by c. Thus 
p is a subspace of g[(f)erQ(iV)) which is invariant under 
the involution A 1— > — A^, where A^ denotes the adjoint 
operator of A with respect to the Hilbert-Schmidt inner 
product on ^zXq{N). Then, q constitutes a Lie subal- 
gebra of 0[(()erg(A^)) which is also invariant under the 
involution A ^ —hJ. By choosing an orthogonal basis 
in \)ZXq(N), this invariance of q translates into a matrix 
representation of g which is stable under X ^ —X'^. 
Then Proposition 1.59 in [HH] implies that g is reduc- 
tive and thus it decomposes into a direct sum of its cen- 
tre 3 and its semi-simple commutator ideal go '■= [fljfl], 
i.e. Q = 3©0o- Since adsu(7v)) is contained in g, the centre 
I is either trivial or R • 1. Thus, similar to Corollary 7.10 
in [HHl, one can show that G (cxpg) is a closed con- 
nected subgroup of GL(f)ero(A^)). Therefore, Theorem 
HLH applies to G. In particular, t and p yield the required 
eigenspace decomposition of Q. Hence we conclude that 
exp(— c) • (expt) = exp(— c) • Kdsu(N) is a Lie subsemi- 



group of GL{\)zxi^{N)) with Lie wedge (— c) © adsu(Ar)- 
Thus the result follows. ■ 

The previous findings suggest the following procedure 
to compute or at least to approximate the Lie wedge of 

(i) Check, whether is self-adjoint (implying pos- 
itive semidefiniteness for F^). This is for exam- 
ple the case, if all Vk in Eqn. lfT3|) are Hermitian 
or, equivalently, if the Kossakowski-Lindblad term 
can be rewritten as a sum of double commutators, 
cf. Eqn. l(22|). 

(ii) If (i) holds, find the smallest cone c containing Tl 
and satisfying the conditions of Theorem IIII.5I 

Note that the above procedure yields but an outer ap- 
proximation of the Lie wedge. In general, further argu- 
ments are necessary to obtain equality. For the generic 
two-level system in [s^l, however, equality can be proven 
as the following result shows. 

Corollary III. 6. Let (S) he a unital Yi- controllable two- 
level system with generic Kossakowski-Lindblad term F^. 
Then, the Lie subsemigroup Ps coincides with 

Ps = exp(-c) • k&su{2) C Co(()ero(2)) 
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(a) (b) 




Figure 2: (Colour online) Steering dynamics of open relaxative systems represented by semigroup actions on a state space 
manifold M: (a) gradient-like method on the reachable set Reach(p) itself; admissible directions are confined to directions 
available in the Lie wedge; (b) optimal control approach as an 'implicit method' on the reachable set Reach(p) brought about 
by a gradient flow on the set of control amplitudes as in Fig. 3 of Ref. [i^. Note that in (b) the entire trajectory at all points in 
time is updated from fe i~> A: + 1 thus exploring more directions than in (a), which may be an advantage over local gradient-like 
methods in open systems. 



where c denotes the convex cone 

c := convjAerLe^ I A > 0, e e Ad5[;(2) } (24) 

contained in the set of all positive semidefinite elements 
in 0[([)ero(2)), cf. Remark \lII.l[ Here, 6^ denotes the 
adjoint operator of 6 with respect to the Hilbert- Schmidt 
inner product on f)ero(2). Moreover, the Lie wedge o/Ps 
is given by (— c) © ad5u(2)- 

Proof. H-controllability of the system implies that 
ad5u(2) is contained in L(Ps). Moreover, for TV = 2 it is 
known that Tl\^ is a positive semidefinite operator 
of g[(f)erQ(2)). Thus Theorem IIII.5I applied to the cone c 
given by Eqn. (|24| yields Ps C exp(— c) ■ Adsu{2)- For 
the converse inclusion, we refer to a standard convexity 
result on Lie saturated systems, cf. [l3|. II 

The geometry of reachability sets under contraction 
semigroups is illustrated and summerised in Fig. [TJ 

In general, it is quite intricate to show that outer ap- 
proximations of the Lie wedge L(Ps) derived from Theo- 
rem lllLSl in fact coincide with L(P5:). To the best of our 
knowledge, no efficient procedure to explicitly determine 
the global Lie wedge of Eqn. (fT4|) does exist. Thus, for op- 
timisation tasks on Reach(p), one currently has to resort 
to standard optimal control methods. A straightforward 
and robust algorithm is mentioned in the final section. 
Moreover, a new approach based on an approximation of 
L(P5:) is sketched. 

IV. RELATION TO OPTIMISATION TASKS 

We follow [Hi in considering optimisation tasks that 
come in two scenarios, see also Fig. [2l (a) abstract opti- 
misation over the reachable set and (b) optimal control 
of a dynamic system specified by its equation of motion 
(e.g. of Kossakowski-Lindblad form). More precisely, an 
abstract optimisation task means the problem of finding 



the global optimum of a given quality function / over the 
reachable set of an initial state p (independently of the 
controls that may drive the system to the desired opti- 
mum). In contrast, a problem is said to be a dynamic 
optimisation task if one is interested in an expHcit (time 
dependent) 'optimal' control that steers the system as 
closely as possible to a desired final state, where 'optimal' 
can be time- or energy-optimal etc. 

In cases where the reachable set Reach(p) can be char- 
acterised conveniently — as, for instance, in closed quan- 
tum systems where it is completely characterised by the 
system Lie algebra so that Reach(p) coincides with the 
system group orbit — numerical methods from non-linear 
optimisation (on manifolds) are appropriate to solve ab- 
stract optimisation tasks on Reach(/9). Details have 
been elaborated in [i^. However, in open quantum sys- 
tems a satisfactory characterisation of the reachable set 
Rcach(/9) — e.g., via Lie algebraic methods — is currently 
an unsolved problem. Thus numerical methods designed 
for optimal control tasks (b) may serve as handy sub- 
stitutes to solve also abstract optimisation tasks (a) on 
Reach(/9) . 

To be more explicit, we consider the Kossakowski- 
Lindblad equation ifTQ]) with controlled Hamiltonian ^ 
in superoperator representation. We are faced with a sys- 
tem taking the form of a standard bilinear control system 
(S) for vec p e reading 

m 

vec p ~ (Ao + ^ UjAj) vecp (25) 

with drift term Aq —i{lN(^Hci—Hj(g)l]y)—TL, control 
directions Aj (E)Hj~Hj (E)tN), and control am- 

plitudes Uj G K, while is given by Eqn. (fTBl) . Then an 
optimal control task boils down to maximising a quality 
functional with respect to some finite dimensional func- 
tion space, e.g., piecewise constant control amplitudes 
(for details see [43| Overview Section). Clearly, one can 



14 



reduce the size of system (|25|) by choosing a coherence- 
vector representation instead of a superoperator repre- 
sentation without changing the principle approach. 

In this context, we would like to point out a remark- 
able interpretation of L(Ps). The method just outlined 
may lead to a (discretised) unconstrained gradient flow 
on some high-dimensional R™. While the 'local' search 
directions (pulled back to state space) are confined to di- 
rections available in the 'local' Lie wedge of Eqn. 
i.e. to the smallest Lie wedge generated by Aq and ujAj, 
Uj G R, the entire method nevertheless allows to vary 
the final point p(T) within an open neighbourhood of 
Reach(p), cf. Fig. [2l[b). In contrast, a gradient-like 
method on the reachable set itself similar to the one for 
closed systems, but with search directions constrained to 
the (local) Lie wedge would in general fail, cf. Fig.[2l[a). 



Outlook: An Algorithm Exploiting the Lie- Wedge 

Yet, combining both methods yields a new approach 
to abstract optimisation tasks: (i) First determine an in- 
ner approximation c of the Lie wedge, (ii) Then, choose 
n G N and define a map from the n-fold cartesian product 
cx---xctoMby(rii,...,f2„) t-^ /(c"" •■•c^^i). Optimise 
this function over the convex set c x • • • x c and increase 
n if necessary. We do expect that the performance of 
such an approach improves the better the approximation 
of the Lie wedge is. In particular, the length of the nec- 
essary products c^^" ■ • ■ will significantly decrease if c 
is a good approximation to L(Ps). Thus even for nu- 
merical aspects knowing the Lie wedge is of considerable 
interest. — With these remarks we will turn to other 
points pertinent in practice. 



Practical Implications for Current Numerical 
Optimal Control 

The above considerations have further implications for 
numerical approaches to optimal control of open sys- 
tems in the sense of the dynamic task (b) of the previ- 
ous section. They provide the framework to understand 
why time-optimal control makes sense in certain WH- 
controllable systems, whereas all other situations ask for 
explicitly taking the Kossakowski-Lindblad master equa- 
tion into account. Consider three scenarios: (i) open 
quantum systems that WH-controllable with almost uni- 
form decay rate, (ii) generic open systems with known 
Markovian (or non-Markovian) relaxation characteris- 
tics, and (iii) open systems with unknown relaxation be- 
haviour. 

In the simple case (i) of a WH-controllable system with 
almost uniform decay rate 7, approximately acts on 
()erQ(Af) as scalar 7I. Now assume that by numerical 
optimal control a build-up top curve giT) (value func- 
tion) of maximum obtainable quality against total dura- 
tion T was calculated for the corresponding closed sys- 



tem with Fi = 0. Moreover, let T* denote the small- 
est time allowing for a quality above a given error- 
correction threshold. Together with the uniform decay 
rate 7 this already provides all information if the quality 
function depends linearly on p{T). Hence determining 
Tl := 'digmax {g{T) ■ e~'^'^} gives the optimal time for 
the desired solution. More coarsely if ~ T*, time- 
optimal controls for the closed system are already a good 
guess for steering a WH-controllable system with almost 
uniform decay rate. 

For case (ii), when the Kossakowski-Lindblad op- 
erator is known, but generically does not commute 
with all Hamiltonian drift and control components, it 
is currently most advantageous to use numerical opti- 
mal control techniques based on the Master equation 
with specific Kossakowski-Lindblad terms as has been 
illustrated in [s^l- The importance of including the 
Kossakowski-Lindblad terms roots in the fact that their 
non-commutative interplay with the Hamiltonian part 
actually introduces new directions in the semigroup dy- 
namics. Likewise, in [IHl, we treated the optimal control 
task of open quantum systems in a non-Markovian case, 
where a qubit interacts in a non-Markovian way with a 
two-level-fluctuator, which in turn is dissipatively cou- 
pled to a bosonic bath in a Markovian way. 

Clearly, the case of entirely unknown relaxation char- 
acteristics (iii), where e.g., model building and system 
identification of the relaxative part is precluded or too 
costly, is least expected to improve by suitable open- 
loop controls, if at all. Yet in [87| we have demonstrated 
that guesses of time-optimal control sequences (again ob- 
tained from the analogous closed system) may — by sheer 
serendipity — be apt to cope with relaxation. In practice, 
this comes at the cost of making sure a sufficiently large 
family of time-optimal controls is ultimately tested in 
the actual experiment for selecting among many optimal- 
control based candidates by trial and error. — Since this 
procedure is clearly highly unsatisfactory from a scientific 
viewpoint, efficient methods of determining pertinent de- 
cay parameters are highly desirable. 



CONCLUSIONS 

Optimising quality functions for open quantum dy- 
namical processes as well as determining steerings in con- 
crete experimental settings that actually achieve these 
optima is tantamount to exploiting and manipulating 
quantum effects in future technology. 

To this end, we have recast the structure of completely 
positive trace-preserving maps describing the time evolu- 
tion of open quantum systems in terms of Lie semigroups. 
On an abstract level, the semigroups of completely pos- 
itive operators may thus be seen as a special instance 
within the more general theory of invariant cones [2i.[89|. 
Here, we have identified the set of Kossakowski-Lindblad 
generators as Lie wedge: the tangent cone at the unity 
of the subsemigroup of all invertible, completely positive. 
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and trace-preserving operators coincides with the set of 
Kossakowski-Lindblad operators. 

In particular, (in the connected component of the 
unity) invertible quantum channels are time dependent 
Markovian, if they belong to the Lie semigroup gener- 
ated by the Lie wedge of all Kossakowski-Lindblad op- 
erators. Moreover, a time dependent Markovian channel 
specialises to a time independent Markovian one, if the 
Lie wedge of an associated semigroup shows the stronger 
structure of a Lie semialgehra. — Likewise, in time de- 
pendently controlled open systems the existence of ef- 
fective Liouvillians that comply with the dynamics given 
by the Master equation is linked to Lie-semialgebra struc- 
tures. 

In view of controlling open quantum systems, reach- 
able sets have been described in the same framework. 
Compared to closed systems, the structure of reachable 
sets of open systems has turned out to be much more del- 
icate. To this end, we have introduced the terms Hamil- 
tonian controllability and weak Hamiltonian controllabil- 
ity replacing the standard notion of controllability, which 
fails in open quantum systems whenever the control re- 
stricts to the Hamiltonian part of the system. For simple 
cases, we have characterised Hamiltonian controllability 
and weak Hamiltonian controllability. These definitions 
also allow for characterising the conditons under which 
time-optimal controls derived for the associated closed 
systems already give good approximations in quantum 
systems that are actually open. In the generic case, how- 
ever, obtaining optimal controls requires numerical tools 
from optimal control theory based on the full knowedge 
of the system's parameters in terms of its Kossakowski- 
Lindblad master equation. 

Finally, we have outlined a new algorithmic approach 
making explicit use of the Lie wedge of the open sys- 
tem. In cases simple enough to allow for a good approx- 
imation of their respective Lie wedges, a target quantum 



map can then be least-squares approximated by a prod- 
uct with comparatively few factors each taking the form 
of an exponential of some Lie- wedge element. 

Since the theory of Lie semigroups has only scarcely 
been used for studying the dynamics of open quantum 
systems, the present work is also meant to structure and 
trigger further developments. E.g., the above considera- 
tions on i-p decompositions may serve as a framework to 
describe the interplay of Hamiltonian coherent evolution 
and relaxative evolution: this interplay gives rise to new 
coherent effects. Some of them relate to well-estabHshed 
observations like, e.g., the Lamb-shift l9C| o r d y namic 
frequency shifts in magnetic resonance |9ll.l92l.l93l|. while 
others form the basis to very recent findings such as 
dephasing-assisted quantum transport in Hght-harvesting 
molecules M, mM,^- 
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